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I. INTRODUCTION 

Non- linear modifications of the Einstein-Hilbert action 
have a long history |l( (for a discussion on recent issues 
see 0, E| and references therein). They have been of in- 
terest for a variety of reasons. For instance, it has been 
claimed, that they could be good renormalizable mod- 
els for quantum gravity 0, [5, 6J. Also some non- linear 
Lagrangians can be chosen with the property that the 
field equations of the metric are second order, these are 
the so-called Lovelock actions 0, SHI El El , which arise 
from dimensional reduction of the Euler characteristic. 
Also string theory predicts an effective gravitational ac- 
tion containing the usual Einstein term plus higher-order 
perturbative corrections in the curvature [T^. Il3| . When 
the Lagrangian is an exact function of the Ricci scalar, a 
mapping exists to the usual Hilbert-Einstein Lagrangian 
with a specific self-interacting real scalar field; the equa- 
tions of motion of the latter being the so called Einstein- 
Klein-Gordon (EKG) system. This idea has been used, 
in different manners to relate well-known scalar inflation- 
ary potentials with pure- gravity higher-order curvature 
scalar Lagrangians El El 111 El El E3 

On the other hand, the EKG system has been studied 
in many situations. In particular, it has been shown that 
there are spherically symmetric solutions that are globally 
regular and asym p toticall y fl at; these solutions are called 
oscillatons El El El El III El E3 El • 

Our aim in this paper is to investigate the mapping of 
oscillatons to non-linear gravity (NLG) theories. We shall 
search for the corresponding objects in these theories, 
and whether we can get interesting information relating 
the results in the NLG theory and those in the EKG 
system. 

We will limit our models to specific scalar potentials 
for which NLG theories can be constructed, and we will 



'Electronic address: octavio@fisica.ugto.mx 
^Electronic address: lurena@fisica.ugto.mx 
t Electronic address: fs@thp.uni-koeln.de 



study the "mapped" objects that may arise in them. It 
is not argued that these gravity models would explain 
all the range of g ravity experiments[27 ^ . Probably, as in 
string theory [13, [13, 28], higher order perturbative cor- 
rections to the Einstein-Hilbert action depending also on 
combinations of the Riemann and Ricci tensors should 
also be present in a more realistic theory, and the mod- 
els considered here would be, possibly, a limit of these 
theories. 

To set the stage for our analysis, we begin by recalling 
some properties of the NLG and EKG systems. For the 
former, we write 



C = y/=~gf{R) 



(1) 



where / is an arbitrary function of the scalar curvature 
R. The field equations derived from are 

f(R)R^-\g^f{R) + {g^n~ V M V„)/'(i?) = 0, (2) 

in which f{R) = j| 0. nnra 

Following a well-known procedure [La Ha . l29| . we now 
establish the connection between this NLG theory and 
the EKG formalism. We first define the new variables 



(3) 



where the conformal transformation between the metrics 
is invertable and then well defined. Inserting this into 
Eq. J2J), we directly obtain the EKG equations 

G>„ := — -g^R = k T^(^) , 



dV_ 



0. 



where 



(4a) 
(4b) 



+ 2V($)], (5) 



is the energy-momentum tensor of the scalar field <& en- 
dowed with a scalar field potential V(<fr). The scalar field 
potential is related to / and R through 

n$) = -^(ii/'-/) . (6) 

2k /' 



2 



Also, K — 8irG (we use units in which H = c = 1, 
and then the Plank mass is mp\ = G 1 / 2 ) and □ = 
g)dfj,(y/— gg^dv) is the covariant d'Alembertian 
operator. Notice that all metric quantities in the EKG 
system will be denoted with a tilde. The corresponding 
Lagrangian is 



R - <f> a <f> , 



2V($) 



(7) 



where we see that the scalar field is minimally coupled 
to the metric tensor g^ v . 

As a final point in this section, we would like to men- 
tion the possible vacuum solutions of Eqs. J2J) known so 
far. It has been shown in|30| that, provided f(R) = 
R + a^R 2 + ■ ■ ■ 44], if a2 > then the only static spher- 
ically symmetric and asymptotically flat solution with 
a regular horizon is the Schwarzschild solution. As it 
can be seen, the statement above is not as strong as the 
Birkhoff theorem that arises in General Relativity, that a 
spherically symmetric gravitational field in empty space 
must be static, with a metric given by the Schwarzschild 
oneHJHi. 

We should emphasize that, as we shall show below, 
the transformation @ will allow us to map known EKG 
oscillaton solutions onto a full NLG theory; that is, the 
transformation will not be performed on the perturbative 
expansion of a NLG theory. 

Although we will exhibit the corresponding perturba- 
tive expansion of the resulting NLG Lagrangian, it is not 
the aim of this work to find solutions to any NLG pertur- 
bative expansion. Nevertheless, the search for oscillaton 
kind of solutions to a NLG theory could be of interest for 
the case of other Lagrangians that, in particular, are ex- 
pressed through a (well defined) perturbative expansion. 

Even though we have shown the formal equivalence 
between two conformally related frames, it should be re- 
minded that this does not imply a physical equivalence 
too. We will make use of the former to find solutions to 
Eq. J2J), and will comment on the latter equivalence in the 
last sections. However, the present manuscript should 
not be seen as a discussion on the trueness of one partic- 
ular frame, but rather as another example that may help 
us to understand some of the (yet hidden) particularities 
of the NLG-frame. 

The sections are organized as follows. In section^ we 
shall present the solutions to the EKG system and many 
of their interesting features. In section IIIII we find the 
corresponding NLG theory and exhibit a way to compare 
the mass-radius relation of the objects formed in both 
systems. Section HVI is devoted to final remarks. 



II. OSCILLATONS: NUMERICAL SOLUTIONS 
OF THE EKG SYSTEM 

Regular and asymptotically flat solutions of the EKG 
equations J3J are called oscillatons, and then we expect 
them to be related to the vacuum solutions of the field 



equations J5J). As an example, we will focus our attention 
in the simplest oscillaton, which arises from a quadratic 
scalar potential of the form V($) = (l/2)m|4> 2 , in the 
spherically symmetric case. Using the polar-areal slicing, 
we write the line element as 



ds 2 = -a 2 (t, f)dP + d 2 {t, f)df 2 + f 2 dQ, , 



(8) 



where a(t, f) and a(i, r) are the metric functions that 
appear on the left hand side of Eq. (|4a|l . 

It is not possible to find analytical solutions to the 
EKG system, but Eqs. Q can be solved numerically in- 
stead. In order to find well behaved solutions, the numer- 
ical solution should be fully time-dependent. The most 
popular solution is due to Seidel & Suen [2TJ, 0j , that 
used Fourier expansions of the functions involved in the 
EKG equations. The method has been refined in |2^.l23| 
to facilitate the numerical solution; and it is the latter 
which is briefly described next to construct the oscilla- 
ton solutions. 

The Fourier series of the fields are of the form 

oo 

V8nG r) = ^^(f)cos(jwi), (9a) 

OO 

a 2 (t,f) = ^a,(f)cos(ju;i), (9b) 



3=0 



a 2 



(t,r) = ^Oj(r)cosO'wi), (9c) 

3=0 



where u) is the fundamental angular frequency of the sys- 
tem. 

Imposing boundary conditions of regularity at the ori- 
gin and of asymptotic flatness, the EKG system becomes 
an eigenvalue problem. For each value of, say, 0i(O), 
there is a set of eigenvalues of w and <Sj(0) for which 
the boundary conditions arc fulfilled (for more details see 
[23|]). The Fourier series have to be truncated by hand, 
and just a few of the Fourier coefficients are taken into 
account. However, all the solutions show convergence: 
the higher the Fourier mode, the less it contributes to 
the series. 

Typical profiles of the Fourier modes of the radial met- 



ric field g r 



are shown in Fig.^ up to j = 12, which 



corresponds to an oscillaton with total mass AI osc — 
0.571(m 2 31 /m 3 >) (here </>i(0) = 2y/2). It should be noticed 
that the time-dependent Fourier terms are confined, and 
that outside a typical radius, say r 9 5 (the radius within 
which 95% of the total mass is contained in), the metric 
matches the Schwarzschild solution, 



a 2 (t, f > r 95 ) 



a- 2 



(t,f> f 9S ) = 1 - 



2GM n 



(10) 



We should add here that the mass of an oscillaton is finite 
due to the fact that the scalar field vanishes exponentially 
for f > f 95 . 

A curious point is that, by construction of the solution, 
$ has only odd multiples of the fundamental angular fre- 
quency lu, and the metric functions a and a have the 
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with respect to R, which yields 
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FIG. 1: The Fourier coeffiecients corresponding to the radial 
metric function g rr = a 2 according to expansion l)9c|l . for an 
oscillaton with total mass M = 0.571(rrip l /m<s>) and for which 
4>i (0) = 2\/2; also shown is its 95% radius rg^ = 9.34771J 1 (For 
oscillatons, m$ becomes the natural unit of distance). The 
Fourier series was truncated at j = 12, and the profiles were 
appropriately scaled to show them all in the same plot. The 
convergence of the solution is manifest. 



even ones. This seems to indicate that the oscillaton in 
Fig. n an d alike are the simplest configurations (the less 
massive) one can construct for the EKG system: the in- 
clusion of all cosine coefficients or of sine terms would 
result in more massive oscillatons |45j. 

An important issue of oscillatons that have not been 
fully tackled is that of their stability. The numerical 
experience so far points out that there exist intrinsi- 
cally stable (S-branch) and unstable (U-branch) oscilla- 
tons [24| ; cf . the more detailed description in catastrophe 
language in [33l l34l |. The U-branch solutions migrate 
away from the equilibrium solution under small pertur- 
bations. Stable oscillatons, on the other hand, do not 
migrate if perturbed a little. Moreover, it seems that 
they play the role of final states an arbitrary scalar field 
configuration evolves towards to. However, there is not a 
definite proof for the stability of oscillatons, but it seems 
that, if oscillatons are not fully stable, they are long lived 
at least [ii|2l Hil . 



III. THE MASS-RADIUS RELATION IN THE 
NLG AND EKG SYSTEMS 



Now, we are to write explicitly the connection between 
the NLG and the EKG systems for the particular example 
of a quadratic oscillaton. Eq. is usually taken as a 
differential equation for f(R) to be solved in terms of the 
curvature R. As can be seen from Eq. J2J), one obtains 
a highly non-linear differential equation, which is very 
difficult to solve. 

This can be ameliorated if we derive Eq. © once again 



2/to/" 



2fV I ■ 1 f> dV R 
^(2/3)k d§ 2k 



(11) 



The obvious and trivial solution, which appears for all 
cases, is f"(R) = 0, which means f(R) = AR+B, with A 
and B arbitrary constants, and also that <E> = const. [ir^. 

It turns out that, in the general case, / and R can 
be given in parametric form in terms of $ by means 
of Eqs. ||BJ| and JIT}. To simplify the calcula- 

tions, we define a new dimensionless variable by ln(x) = 
V(2/3)k $. Thus, 



R 



2kqx 2 
Akqx 



V{x) + x 
V(x 



dV(x) 



dx 
x dV{x) 
2 dx 



(12a) 
(12b) 



In this form, it is easily seen that the Liouville the- 
ory, for which V($>) = Vbe A *, where Vb and A are con- 
stants, is one of the scalar-tensor cases that can be solved 
exactly [II EEI. However, it seems that the Liouville 
theory does not provide regular nor asymptotically flat 
solutions of the EKG system, see for instance^ |47|. 

On the other hand, for the quadratic case V($) = 
(l/2)m|$ 2 , Eqs. O ar e rewritten as 



/ = -m|x 2 ln(a;)[ln( a ;) + 2] , 



(13a) 

R = 3m|iln(a;) [ln(x) + 1} . (13b) 
A Taylor expansion of Eqs. (|P3*|) around /(0) = (x = 1) 



gives 



f(R) = R 



1 



6m| 



R 2 



1 



18m; 



-i? 3 + 0(R 4 



(14) 



Hence, a quadratic oscillaton corresponds to an i?-regular 
NLG theory. 

It has been argued before that Lagrangians of the 
form (|14fl accomplish some desirable features (3^. For in- 
stance, the presence of an i?-term ensures that the EKG 
exists near Minkowski spacetime, and the existence of 
an a2i? 2 -term ensures regularity of the conformal trans- 
formation to flat space; moreover, if 02 > then the 
Minkowski space is a stable ground state solution. It 
also seems that curvature corrections at all orders are 
essential in order to regulate gravity [2^. 

It should be remarked that NLG theories with char- 
acteristics as mentioned in the above paragraph corre- 
spond to a very particular class of scalar potentials and 
probably are a limit of a more general theory From 
Eqs. I|13|) . it can be seen that Lagrangians of the form 114|) 
correspond to the case in which V(0) — ^g-(O) = 0, i.e., 
the scalar potential V($>) has a critical point at <f> = 0. 
Explicitly, the first Taylor coefficients for this type of po- 
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tentials, calculated from Eqs. 112(1. are 
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= X 
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18 [V"(0)] 



where ' = -4- and ' = -A . Higher Fourier terms depend on 
higher derivatives of the scalar potential in a complicated 
way, so that we do not show them here. However, it is 
easily seen that 02 > only if the scalar field is also 
massive, i.e., m\ = V"(0) > 0. 

What else can we say of the NLG theory in 1(13|) ? First, 
the line clement in the NLG theory is spherically sym- 
metric and fully time-dependent, properties preserved by 
the conformal transformation • The metric fields have 
Fourier expansions of the form 



a 2 = e -V( 2 /3)-o* 5 2 = J2 aj (f) COs(jLjt) , (15S 
3=0 



a = e 



V(2/3)«o* a 2 = aj (f) COB(jWt) , (15b) 
3=0 



where, in contrast to Eqs. 10 of an EKG-oscillaton, all 
the Fourier coefficients, odd and even, are non-trivial. 
Here, it is the scalar field <& which provides the odd co- 
efficients. 

Second, the exterior form for both metrics g and g 
is the Schwarzschild solution Q1U[I. In consequence, the 
ADM mass is the same in both the NLG and the EKG. 

Therefore, this theory permits the existence of regu- 
lar objects made of pure gravity, which are, in addition, 
vacuum solutions. Hence, there may be in some cases 
a com plem entary second part of the theorem mentioned 
before 30] . We have shown here that there is a non-linear 
Lagrangian of the form f{R) = R + a.2-R 2 + • • • , with 
&2 > 0, that has two type of asymptotically fiat vacuum 
solutions: Schwarzschild (static), and that corresponding 
to oscillatons (time-dependent). It would be interesting 
to investigate whether the above statement is true for 
any massive real scalar field, since it has not been proved 
completely that the EKG equations have regular stable 
solutions only in the case m% > 0; though this seems to 
be the case, see for instance [37j. 

We observe that the form of all of the coefficients in the 
Taylor expansion of the exact result l|13(l is determinant 
for the existence of oscillaton solutions; had we started 
with the approximate expansion i|14|) . the results could 
have in general differed from the results obtained in here, 
even approximately. 

On the other hand, it should be noticed that the sim- 
plest solution to Eqs. J2J) is the metric g^ whose coeffi- 
cientes are expanded in a Fourier series of the form (1151) . 
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FIG. 2: (Top) Comparison of the radial metric function ob- 
tained from the EKG and the NLG cases for the same case 
shown in Fig. Q It can be seen that the NLG solution is 
asymptotically flat and coincides with the EKG solution for 
r — » 00; but it is conformally flat at f = 0. (Bottom) M 
vs f max plots for the EKG and NLG solutions, the configu- 
ration corresponding to the leftmost (rightmost) point of the 
plot has ^/7^$(0,0) = 0.8 (0.01). The NLG-oscillatons are 
less compact that their EKG-counterparts, so that they are 
distinguishable in principle, see text for details. 



As we said before, this metric is spherically symmetric 
but is not written in its so-called standard form. An ob- 
server in the NLG frame would, however, notice that the 
coordinates i, f are its usual Schwarzschild coordinates 
far from f = 0. This would indicate her/him that there 
is an exterior (static) vacuum solution and an interior 
(time-dependent) one. 

This observer would also find more handy to continue 
using (i, r) as her/his metric coordinates, and then to 
preserve g^ v in its simplest form. Thinking of this pos- 
sibility, we show a comparison between the radial metric 
functions in the EKG and NLG cases in Fig. [21 We also 
plot the corresponding M vs f max graphs for both cases, 
where f max is the position of the maximum of the radial 
metric function (|15b|) at t = 0. There is a maximum mass 
for the NLG case, which could be interpreted as an indi- 
cation of the existence of Stable and Unstable branches, 
as in the case of the EKG theory. 



5 



IV. FINAL REMARKS 

If the NLG theory in (|13[) were the fundamental gravity 
theory, there would be many dark pure-gravity objects, 
corresponding to oscillatons, in the universe. For a given 
mass, the object in the NLG theory is bigger than its 
EKG counter part (Fig.[2J), so that one could, in principle, 
discriminate between them. Ultimately, by these means 
we could determine which system, the NLG or the EKG, 
is being measured. 

On the other hand, we could ask at this point: how 
can these objects be formed from the NLG point of view? 
As we mentioned before, in the EKG, a self-gravitating 
scalar field evolves and settle-down onto an S-oscillaton. 
Because of the spherical symmetry of the system, the 
only mechanism of relaxation permitted is the emission 
of scalar radiation (scalar field matter). This pro cess has 
been dubbed gravitational coolina^lL l24l l25l l38||. which 
has been shown to be a very efficient mechanism for the 
relaxation of self-gravitating scalar fields. 

On the other hand, any initial scalar field configuration 
in the EKG frame can be conformally transformed into 
a pure gravity configuration in the NLG. Moreover, also 
the complete evolution of the EKG system can be con- 
formally transformed and followed in the NLG. As the 
system is spherically symmetric also in the NLG, gravi- 
tational radiation is forbidden, and so there should exist 
a mechanism for the relaxation which would be the "con- 



formal" partner of the gravitational cooling. 

The only explanation we foresee is that, in the higher 
derivative NLG frame, gravity must be allowed to have a 
spin-0 component, which would provide the channel for 
gravitational cooling. In fact, the existence of scalar grav- 
itational waves has been studied before, together with the 
idea that actual interferometers built for the detection of 
gravitational waves can also detect a scalar component 
of gravitational radiation[3!j ^(j ■ 

It would be interesting to determine the scalar grav- 
itational waves emmited for the system discussed here, 
for which one would expect that the non-linear-terms 
would, in a certain manner, take the role of the scalar 
field. The procedure would be the full numerical evolu- 
tion of Eqs. J2J). For this, it would also help the formal 
equivalence between the NLG and the EKG, as the nu- 
merical evolution of the metric fields a and a can be 
mapped directly from that of a, a and These calcula- 
tions are in progress and will be reported elsewhere. 
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